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ABSTRACT 

We estimate the mean dust-to-gas ratio of MgH absorbers as a function of rest equiv- 
alent width Wq and redshift over the range 0.5 < z < 1.4. Using the expanded 
SDSS/HST sample of low-redshift Lyman-a absorbers we first show the existence of a 
8a correlation between the mean hydrogen column density (iVm) and Wo, an indica- 
tor of gas velocity dispersion. By combining these results with recent dust-reddening 
measurements we show that the mean dust-to-gas ratio of MgH absorbers does not 
appreciably depend on rest equivalent width. Assuming that, on average, dust-to-gas 
ratio is proportional to metallicity, we find its redshift evolution to be consistent with 
that of L* galaxies from z — 0.5 to 1.4 and we show that our constraints disfavor dwarf 
galaxies as the origin of such absorbers. We discuss other scenarii and favor galactic 
outflows from ~ L* galaxies as the origin of the majority of strong MgH absorbers. 
Finally, we show that, once evolutionary effects are taken into account, the Bohlin 
et al. relation between Ay and Njj is also satisfied by strong MgH systems down to 
lower column densities than those probed in our Galaxy. 
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1 INTRODUCTION 

A successful theory of galaxy formation must not only ex- 
plain the properties of the luminous parts of galaxies but 
it is equally important to account for the material seen in 
absorption against background sources. However, while the 
connection between stro ng absorbers and galaxies was re- 
alized a long time ago (|Bahcall fc Spitzerl Il969i ; iBergeronl 
1986), our physical understanding of absorption selected sys- 
tems has not yet reached the maturity of current models 
addressing the emission properties of galaxies. In partic- 
ular, the nature of the structures probed by some of the 
strongest absorption features is n ot yet understood (e.g. 
IWolfe. Gawiser. fc Prochaskall2005l ). 

Exploring the existence of scaling relations for absorbers 
may provide us with useful insight to constrain the physical 
conditions and establish theoretical models for these systems 
and their associated galaxies. In this paper we investigate 
the properties of the dust-to-gas ratio of MgH absorbers, 
i.e. the strongest metal line detectable in optical spectra at 
z < 2 and in some c ases a tracer of Damped Lyman-a ab - 
sorbers (DLAs) (e.g. IChurchill. Kacprzak. fc Steidelll2005l ). 
We first present new correlations involving gas velocity dis- 
persion, hydrogen and dust column densities, and show how 
the knowledge of the dust-to-gas ratio sheds light on the 
nature of strong MgH absorber systems. 

The dust-to-gas ratio is one of the basic properties 



of the ISM and IGM but our knowledge of this quantity 
is largely based on our own Galaxy, some of its satellites 
and a few objects with z > using f or example mul- 
tiple images in strong ly lensed systems l|Zuo et al.l 1 19971 : 
iDai fc Koc hanek 2008) or absorbers in front of gamma ray 



bursts (|Ellison et al.l l200fj ). Interestingly, the Milky Way 



presents a characteristic value of the extinction per H atom. 
Using hydrogen Lyman-a and H2 absorption lines to de- 
termine the total H column de nsities along star sight lines, 
iBohlin. Savage, fc Drake] (|l978l ) found that the extinction in 
the U-band follows 



0.53 



10 21 cm" 2 



mag 



for Rv = 3.1 , 



(1) 



with a scatter about the mean of about 30% and where 
Rv = Av/E(B — V). This relation was initially observed 
over the range of column densities 10 20 < N H < 3 x 10 21 
cm -2 and then extended and confirmed up to Nh ~ 5 x 10 22 
cm" 2 where molecular hydrogen plays an important role 
|Snow. Rachford. fc Figosklll2002l ). The linearity and slope 
of the above relation carry information on the complex 
mechanisms responsible for the formation and destruction 
of dust grains. It is interesting to investigate whether such 
a relation holds at lower column densities (where ionization 
corrections could be important) and/or in different environ- 
ments such as larger galactic radii and higher redshifts. 

By extrapolating the above relation, column densi- 
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Vladilo. Prochaska fc Wolfe! (I2008T) and Mgll absorbers 
York et all (|2006h . I Wild et al.1 (|2006l ) and IMenard et all 



ties with Nh < 10 20 cm -2 are expected to be associ- 
ated with E(B-V) values smaller than a fewxlO -2 mag- 
nitude. Measuring such an effect along individual star sight 
lines becomes challenging and a statistical approach is re- 
qu ired. Such a technique was recen tly used with DLAs 
by 

bvJLf 

(2008). In particular, the latter authors constrained the 
reddening induced by Mgll absorbers down to E(B-V)~ 
5 x 10~ 3 mag by analyzing ~ 7000 intervening systems. 

By combining reddening constraints of Mgll absorbers 
with measurements of hydrogen column densities, the 
present analysis will allow us to obtain constraints on dust 
properties in new regimes, i.e. at z ~ 1 and on large scales 
(~ 50 kpc) around galaxies. In particular we will show that 
the in ferred dust-to-gas ratio of Mgll i s similar to that given 
by the lBohlin. Savage, fc Drake! (1 19781 ) relation extrapolated 
down to lower hydrogen column densities. 

Throughout the paper we will simply refer the dust-to- 
neutral gas ratio as dust-to-gas ratio unless stated otherwise. 
The amount of dust will be quantified by either E(B-V) or 
Ay. The Mgll rest equivalent width, Wo, denotes the 2796 A 
transition. 



2 THE DATA 

The data used in this analysis come from two sources: (i) 
"the expanded SDSS/HST sample of low-redshift Lyman-a 
absorbers" compiled by Rao et al. (2006) and (ii) the 
characterization of the re ddening effects induc ed by Mgll 
absorbers in the SDSS by IMenard et all (|2008D . We briefly 
describe them below. 

The expanded SDSS/HST sample of low-redshift 
Lyman-a absorbers consists of 197 systems and represents 
the largest sample of UV-detected DLAs ever assembled. It 
has been compiled from UV spectra of QSOs with Mgll sys- 
tems optically identified in the range 0.11 < z < 1.65. The 
Mgll absorber systems were selected from various sources 
in the literature (See Rao. et al 2006 for references). Most 
of the UV data were obtained with observing programs on 
HST led by S. Rao as well as archival data. 

In Fig. Q] we show the distribution of neutral hydrogen 
column density, Nhi, as a function of Mgll rest equivalent 
width, with 0.45 < Wo < 3.3 A. The arrows show upper 
limits on Nhi- They represent about 3% of the data points 
and can safely be neglected in the present analysis. An im- 
portant point to emphasize is the large scatter in Nhi (up 
to three orders of magnitude) even at a fixed value of Wo- 
The color-coding of the points refer to the measured rest 
equivalent width of Fell at 2600 A. Empty circles are used 
for upper limites. 

The second set of observational results used in our anal- 
ysis comes from the recent dust-reddening analysis done by 
IMenard et"aH (|2008l ). Using almost 7,000 Mgll absorbers de- 
tected in SDSS quasar spectra these authors have measured 
the mean color excess E(B-V) induced by these systems, as 
a function of rest equivalent width and redshift. Such con- 
straints will be used to compute the dust-to-gas ratio of Mgll 
absorbers. 

It is useful to recall that the majority of Mgll absorp- 
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Figure 1. Distribution of hydrogen column densities Njji of 
Mgll selected absorbers, as a function of rest equivalent width 
W 2796 . The data originate from the expanded SDSS/HST sam- 
ple of low-redshift Lyman-a absorbers compiled by Rao et al. 
(2006). Typical errors are about 0.1 dex in Nhi and about 0.07 A 
in Wq 796 . The arrows show upper limits in Nhi , the color of the 
points refer to the measured rest equivalent width of Fell, Wq® 00 . 
Empty circles denote upper limits of Wq 600 . The green, red and 
orange points with error bars represent the arithmetic mean, ge- 
ometric mean and median in each bin, respectively. 



tion lines with Wo > 1 A are saturated (|Nestor et al.ll2005T l 
and that, therefore, the rest equivalent width of the lines 
provides us with an estimate of the gas velocity dispersion. 
For a completely saturated line, 1 A corresponds to Av ~ 
107 km/s. Using high-resolution spectroscopy the absorption 
is often seen to originate from several velocity components 
implying that this value is only a lower limit on the total 
velocity dispersion of the system. E mpirically it h as been 
observed that Av ~ 120 km s _1 A" 1 (|Ellisonll2006l . Fig. 3). 
This value will be used below for physical interpretations. 



3 THE DISTRIBUTION OF HYDROGEN 
COLUMN DENSITIES 

The expanded SDSS/HS T sample o f low- redshift Lyman- 
a absorbers was used bv lRao et al.l (120061 ) to quantify the 
amount of neutral hydrogen probed by Mgll absorbers at 
z < 1.65 and to estimate the mean hydrogen density Qh ■ To 
do so, these authors computed the mean HI column density 
as a function of Wo (Mgll): 
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i N 

(N m )(W ) = -J2 Nl »>*> ( 2 ) 

and found (N H i> ~ 10 21 cm -2 for the range Wo > 0.6 A. 
In Figure [1] we show the arithmetic mean of Nhi in green, 
using a logarithmic binning in Wo- We have used 500 boot- 
strap samples in order to estimate the errors on the mean in 
each bin. We recover the lack of corr elation between mean 
hydrogen column density reported by [Rao et all (I2006D . For 
Wo > 0.6 A we find (Nhi) ~ 3 x 10 20 cm" 2 . 

While the use of the arithmetic mean is needed to quan- 
tify the amount of neutral hydrogen probed by Mgll ab- 
sorbers, which is important to constrain the evolution of HI 
through cosmic time, it ends up being sensitive to a small 
fraction of the data points, having the largest Nh values. 
Indeed, the distribution P(Nhi\Wo) is highly asymmetric 
and spans about three orders of magnitude in the present 
dataset. The value of (Nhi) is therefore driven by a small 
fraction of the data points. Even if it is computed from a 
sample of ~200 objects, it might effectively suffer from small 
number statistics. □ 

More important, the arithmetic mean of Nhi does not 
provide relevant information regarding the majority of Mgll 
absorbers and is therefore not suited to extract underlying 
correlations between the parameters Nhi and Wo. However, 
additional information can be extracted and we first con- 
sider the geometric mean of Nhi'- 

(Nhi) b (W ) = 10< logJV *">. (3) 

We compute it using the same binning in Wo and the errors 
are again estimated with bootstrap resampling. The results 
are shown in Figure[T]with the red data points. We can now 
detect a strong correlation between the hydrogen column 
density Nhi and the Mgll rest equivalent width Wo and a 
simple power-law fit gives 

(Nm) (Wo) = C g {WoT 3 (4) 

where the subscript g denotes a geometric mean, 
C g = (3.06 ± 0.55) x 10 19 cm -2 and a g = 1.73 ± 0.26. This 
correlation holds over one order of magnitude in Wo- Its 
significance, quantified by that of the slope, is greater than 
6cr. 

To demonstrate the robustness of the above result, we 
repeat the procedure and estimate the median of Nhi as 
a function of Wo- The results are shown in Figure [1] with 



Table 1. Scaling parameters 
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orange points. They show consistency with the results ob- 
tained using the geometric mean. Similarly, a power-law fit 
gives 

med [Nhi] (W ) = C m {W ) am (5) 

with Cm = (2.45 ±0.38) X 10 19 cm -2 and a m = 2.08 ±0.24. 
The correlation is now detected at 8.7o\ Consistent fitting 
parameters are obtained using both estimators which show 
that, even if there is a large scatter between Nhi and 
Wo, there exists an underlying and well defined correlation 
between them. Both the geometric mean and the median 
allow us to measure a signal coming from the majority of 
the systems, and not from a few outliers. It is interesting 
to mention that the correlation weakens as we increase the 
lower limit for Nh, but it is not severely affected if we 
decrease the higher limit of Nh (the opposite statement 
would apply to the arithmetic mean). Our results indicate 
that, typically, Nhi is roughly proportional to the square 
of W . 

Similarly, we can quantify the relations between Nhi 
and Fell rest equivalent width. By doing so we should keep in 
mind that the sample used in this analysis is Mgll-selected. 
The corresponding Fell-related correlations are valid in this 
context only. 

We have applied the three estimators introduced previ- 
ously and summarize the results in Table [2] The correlation 
between the two parameters is detected at a high signifi- 
cance (8 and 12<r) and is closer to a simple proportionality Q 
. Such a behavior is expected as the 2600 A Fell absorption 
line has a lower oscillator strength that the Mgll doublet. It 
is therefore less subject to saturation and carries more infor- 
mation about the column density of the absorbing system. 



1 The same problem affects estimates of the cosmological density 
of neutral gas, Q g (z), as its estimation involves a similar average 
(Lanzetta et al. 1991): 

H fim H Ei N H i,i 
il g( z l = — ■ 



Pcrit 



AX 



where fi is the mean molecular weight of the gas, vtih is the 
mass of the hydrogen atom, p cr u is the current critical mass 
density and AX is the absorption distance path. This effective 
small number statistic might be at the origin of the discrep- 
ancies in Qhi currently debated in the literature. Based on a 
sample of ~ 2 Mg ll selected systems with Nhi measurements, 
iPeroux et al.l ||2004T ) already reported such an effect. 



4 DUST-TO-GAS RATIO 

The presence of dust associated wit h Mgll absorbers 
has been reported by several authors dMenard Sz Perouxl 
20031; IWang et alj|2004 iKhare et al.ll2005l;lYork et al.ll2006l ; 
Wild et al.ll2006l ). Recently. iMenard et alj i|200Sf ) analyzed 
close to 7,000 strong Mgll absorbers and quantified the 



2 We also note that a double power-law is a more accurate repre- 
sentation of the data in the (Nhi, Wo(FeI I)) plane. This is cer- 
tainly an interesting feature to address but is beyond the scope 
of this paper, focusing on Mgll absorption. 
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amount of reddening as a function of rest equivalent width 
and redshift. In particular, they find a scaling relation be- 
tween the amount of reddening and the Mgll rest equivalent 
width: 







(6) 



(E(B-V) lest )(W ,z)=Cx (j^-j (1 + z) 

where a = 1.88±0.17, P = -1.1±0.4 and C = (0.60±0.07) x 
10~ 2 mag. This relation has been constrained over the range 
1 < Wo < 6 A. The similarity of Eq. |S]and Eq. 2] is striking 
and suggests that the average dust-to-gas ratio of Mgll se- 
lected systems does not strongly depend on rest equivalent 
width. We now quantify this statement by measuring the 
dust-to-gas ratio in the range of Mgll rest equivalent widths 
common to both datasets, i.e. 1 < Wo < 3.3 A. 

To do so we h ave rerun the reddening analysis by 
I Menard et al.l (|2008h selecting only Mgll absorbers with 
0.4 < z < 1.4 so that the redshift distributions of the two 
datasets are similar. It is important to note that, the me- 
dian and arithmetic means give consistent results for the 
reddening analysis 0. This is probably due to the fact that 
high dust column density systems cannot be detected in the 
SDSS survey due to the extinction bias and therefore the tail 
of the redden ing distribution cann ot be probed. Richards et 
al. (2003) and lMenard et all i|200Sl ) showed that quasars red- 
dened by a color excess E(B — V) greater than about 0.3 
magnitude can hardly be selected. Fortunately, this effect is 
not expected to str ongly bias the reddening constraints as 
lEllison et all (|2004 ) showed, using radio-selected quasars, 
that no more than 20% of optically-selected quasars are 
missed due to extinction effects associated with strong Mgll 
absorbers. 
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Figure 2. Upper panel: geometric mean and median estimates of 
Njjj as a function of Mgll rest equivalent width (upper x-axis) 
or gas velocity width (lower axis). Middle panel: Mean redden- 
ing induced by Mgll absorbers from lMenard et all j200Sh . Lower 
panel: The dust-to-gas ratio Ay /Nhi which does not significantly 
depend on the Mgll gas velocity width. 



4.1 Dependence on Wo 

We present our measurements of the dust-to-gas ratio of 
Mgll absorbers in Fig. [2] The upper panel show the variation 
of the geometric mean and median of Nhi as a function 
of Mgll rest equivalent width (upper axis) or Av(Mgll), 
the gas velocity dispersion (lower axis). The lines show the 
fitted power-laws described in Table 1. The middle panel 
shows the observed reddening values E(B — V) for Mgll 
absorbers selected with 0.4 < z < 1 .4. The solid line is 
the fitting formula (Eq.|6j) proposed bv lMenard et al.l (120081 ) 
and simply evaluted at the mean redshift of the sample. To 
present the results in convenient units we convert E(B — V) 
reddening values into visual extinction Ay using a small 
Magellanic cl oud (SMC) e xtinc tion c urve, i.e. Ry = 3.1 a s 
motivated by lYork et all (|2006l ) and IMenard et al l i|200Sl ). 
We use the ratio between < Ay > and < N(HI) > (directly 
obtained from the data) as an estimate of the mean dust-to- 
gas ratio of Mgll absorbers Q. For systems with 1 < Wo < 



3 The geometric mean cannot be computed as the signal origi- 
nates from positive and negative fluctuations. 

4 Ideally, the average dust-to-gas ratio < Ay /N(HI) > can be 
estimated by weighting the individual reddening measurem ents 
by 1/N(HI), as done by I Vladilo. Prochaska. fc Wolf j ll2008h . In 
the present study, the number of objects for which N(HI) is avail- 
able is unfortunately too small to provide us with a reddening 
detection. 



3.3 A, we find using the geometric mean or the median of 
N(HI): 



{Ax 



(N (HI)) 



= 3.0 ± 0.6 x 10~ 22 mag cm 2 



(7) 



at a mean redshift close to unity. This value is less than 
a factor two lower than that of the Milky Way. It is inter- 
esting to note that while the shape of the mean extinction 
curve of Mgll absorbers is consistent with that of the SMC , 
i.e. does not presen t the 0.2 urn bump (|York et al 1 120061 : 
IMenard et ail liboS). the mean dust-to-gas ratio is substan- 
tially higher than that of the SMC. 

Interestingly, the mean dust-to-gas ratio does not ap- 
pear to be a strong function of Wo- While both the mean 
hydrogen and dust column densities vary by more than an 
order of magnitude over this range, the mean dust-to-gas 
ratio is consistent with being constant and is found to vary 
by less than a factor ~ 2 in this interval for the current 
sample. In comparison, the dust-to-gas ratios of the Milky 
Way and the SMC differ by a factor ~ 8.5. A power-law fit 
to the observed dust-to-gas ratio as a function of Mgll rest 
equivalent width gives 



(iV(HI)} 



oc (Wo) 



79 = 

7™ 



0.5 ±0.7 
= 0.4 ±0.8 



(8) 



A large scatter in the relation between dust-to-gas ratio and 
Mgll rest equivalent width may exist if we consider individ- 



© 0000 RAS, MNRAS 000, 000-000 



HI content and dust-to-gas ratio of Mgll absorbers 5 



ual systems, however our results show that the mean dust- 
to-gas ratio of Mgll-selected systems is not a strong function 
of gas velocity dispersion. The lack of strong correlation is 
an interesting property. 

If we assume that the dust-to-gas ratio is, on aver- 
age, proportional to metallicity, the trend reported in Eq. 
[8] indicates that the mean metallicity of Mgll absorbers 
is a weak function of rest equivalent width. The linear- 
ity between dust-to-gas ratio and metallicity is expected 
as dust is formed from metals and the mass of metals is 
equal to the product of the metallicity and the gas mass. 
This linear correlation has been shown in nearby galaxies 
|lssa. MacLaren. fc Wolfendald Il99d; IrBoissier et alj 120041 ) 
or m high redshift DLAs (jVladilo et all 120061 ) over several 
orders of magnitude in column densities. 

Given the existen ce of a stellar mass-metallicity relation 
jTremonti et al.ll2004l ). the apparent lack of correlation be- 
tween the mean dust-to-gas ratio and the velocity dispersion 
of the gas implies that the mean mass of galaxies giving rise 
to strong Mgll absorption does not play a dominant role 
in determining Wo- In other words, the velocity dispersion 
of the gas does not reflect the gravitational potential of 
the system. By analyzing the correlation between the mean 
luminosity of Mgll absorbing galaxies and Wo, several 
authors have reported similar results. Using a sample o f 
58 Mgll absorber-galaxy associations, S teidel et al.1 (|l997l ) 
did not find any correlation between galaxy luminosity 
and absorber rest equivalent wid th. Similar results were 
also found by IZibetti et al.1 ([2007) who stacked the images 
of ~ 2500 SPS S quas ars with strong Mgll absorbers and 
iKacprzak et al.l (|2007l ) who reported a lack of correlation 
between Mgll absorption and galaxy morphology from HST 
data. 

The present analysis allows us to probe only a factor 
~ 3.3 in Mgll rest equivalent width, which corresponds to 
velocity widths in the range 120 < Av < 400 km s _1 . The 
number of absorbers per unit rest equivalent width and red- 
shift is given by Nestor et al. (2005): 



-21 



N* w o 
dN/OWo = —e~ — 



(9) 



with the maximum likelihood values W* = 0.702 ± 0.017 A 
and N* = 1.187 ± 0.052. The steep exponential decline im- 
plies that the incidence of systems with Wo — 1 A is about 
30 times higher than that of systems with Wo — 3.3 A. 
The sharp decrease of dN/dW-o may be due to a transient 
nature of the absorbing gas, such as outflows triggered by 
star formation. 

These results may, at first, appear to be in contrast with 
direct measurements o f metallicities and M gll rest equiva- 
lent widths reported bv lMurphv et al' (2007). These authors 
found (Z) ~ (1.69 ± 0.20) log(Wn) + est which indicates 
a substantially steeper dependence on Wo- However, their 
sample selection significantly differs from ours: their study is 
not based on Mgll selected absorbers in general, but focuses 
only on systems with both detections of hydrogen column 
density and a metallic absorption line from a volatile ele- 
ment. As a result of this additional selection criterion, the 
distribution of points in the {Nhi ,Wo) plane differs from 
the generic one compiled by Rao et al. (2006) and presented 
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Figure 3. The mean dust-to-gas ratio of strong Mgll absorbers 
as a function of redshift. For comparison we show an estimate 
of the expected dust-to-gas ratio of ~ L* galaxies from Dave & 
Oppenheimer (2007) as well as the valu es of the Milky Way, the 
LMC and the SMC l|Gordon et al.ll2003h . 



in Fig. \T\ For example, at Wo < 1 A the region of the plane 
with Nhi < 10 21 cm -2 is almost unpopulated in their case 
while it is where most of the points lie in the Rao et al. sam- 
ple. Hence, the metallicity trend reported by Murphy et al. 
only applies to a specific sub-population of Mgll absorbers. 

We postulate that, similarly to the shallow dependence 
found between dust-to-gas ratio and Wo, a weak correla- 
tion is expected between metallicity and Wo providing that 
absorbers are selected only on their Mgll rest equivalent 
width. Such a meas urement can be performed by using com- 
posite spectra (e .g. iNestor et al.ll2003l ; iTurnshek et al.|[2005l : 
I York et al.ll2006T l which allows one to detect weaker absorp- 
tion lines and estimate mean metallicities. 



4.2 Redshift evolution and link to L* galaxies 

Having shown that the dust-to-gas ratio of Mgll absorbers 
does not strongly depend on rest equivalent width (over the 
range 1 < W < 3.3 A), we now investigate its evolution as a 
function of redshift. Following the above procedure, we have 
divided both datasets into two redshift bins with (zi) ~ 0.7 
and {22) — 1.2 and measured the mean Av /Nhi{z) as de- 
fined above. We present the results in Fig.[3]and the numer- 
ical values are given in table 2. As can be seen, the observed 
dust-to-gas ratio of strong Mgll systems is, on average, sim- 
ilar to that of the Milky Way and significantly higher than 
that of the SMC, even at the highest redshift we can probe. 
The data also suggest a trend of decreasing dust-to-gas ratio 
with increasing redshift. 

As done previously, assuming that the dust-to- 
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Table 2. Dust-to-gas ratio of Mgll absorbers with 1 < Wq < 3 A 



Av/N ffl [lO" 22 mag/cm 2 ] 


(z) = 0.7 


4.20 ±0.14 


(z) = 1.2 


3.10 ±0.74 



gas ratio is, on average, proportional to the metal- 
licity (see section I4.1|l . we can attempt to compare 
the observed trend to available models. To do so, we 
first consider that a solar metallicity corresponds to a 
Milky Way dust-to-gas ratio dBohlin. Savage, fc Drake! 19781 ; 
llssa, MacLaren. fc Wolfendalelll99cK ~ 



log 



A v , . 

IN HI 



1os( Nh7 



(10) 



The above estimate strongly depends on the averag- 
ing procedure used to define the mean metallicity. Mod- 
elled evolutions of ([ Z/H1) (z) have been explored by 
iDave fc Qppenheimerl (|2007l ). These authors investigated 
the cosmic metal budget in various phases of baryons es- 
timated from cosmological hydrodynamic simulations which 
included constrained models for enriched galactic outflows. 
The solid line in Fig. [3] shows their estimation of the varia- 
tion of mean metallicity as a function of redshift for a star- 
formation weighted estimator, expected to be a representa- 
tive value for galaxies selected in emission. We can see that, 
over the entire redshift range available, the dust-to-gas ra- 
tio of Mgll-selected systems is in agreement with the model 
of the metallicity evolution using a star formation-weighted 
estimator (representative of L* galaxies) but not with the 
N(HI)-weighted one (corresponding to DLAs). Even at the 
highest redshifts probed by our dataset, the mean dust-to- 
gas ratio of Mgll absorbers is significantly higher than that 
of the SMC at z — 0. For comparison, we also show the ex- 
pected evolution of an JV^i-weighted estimator of ([Z/H])(z) 
with the dashed line. As can be seen, such a trend is not in 
good agreement with the data points. 

These results show that the mean dust-to-gas ratio of 
strong Mgll absorbers is consistent with that of L* galax- 
ies but not with that of substantially smaller systems, such 
as the SMC. Given that metallicity is expected to decrease 
with increasing redshift, an LMC-type dust-to-gas ratio is 
also disfavored. It suggests that, on average, strong Mgll 
systems are associated with ~ L* galaxies (which does not 
prevent the existence of a large scatter around this relation 
if individual systems are considered). 

Associations between strong Mgll absorbers and ~ L* 
galaxies have already been made obse rvationally. A s men - 
tioned above. Isteidel et all (|l997n and IZibetti et all (|2007n 
showed that the mean luminosity of Mgll absorbing galax- 
ies is about 0.8 L* . However, while such associations clearly 
gathered information on the link between strong Mgll ab- 
sorbers and ~ L* galaxies, they have not been able to pro- 
vide strong constraints on the nature of these systems. Below 
we show how the knowledge of the mean dust-to-gas ratio 
may shed light on the origin of these systems. 



4.3 Nature of the absorbing gas 

The origin of Mgll absorbers has been a matter of de- 
bate since their discover y. Various scenarii have bee n pro- 
posed: infehMiigmaterial l|Mo fc Miralda-Escudd 19961). out- 
flows l|Bond et alfeooj ) , orbiting dwarf galaxies ( York et all 
1986), etc. Certain studies have proposed a relation between 
strong Mgll absorbers and star forma tion: from the observed 
redshift distribution of dN/dz (e.g. Guillemin fc Bergeron! 



1 19971 : iProchter. Prochaska. fc Burled 120061 ) or from the ex- 
istence of a correlation between Mgll rest equivalent width 
and star formation rates estimated from broad band colors 
IZibetti et ail (2007). While some indication of the outflow 
phenomenon exists in certain cases, two different origins re- 
main: outflows can equally well originate from the bright 
galaxy usually found at ~ 50 — 100 kpc from the absorber 
or from the satellite galaxies orbiting the parent dark matter 
halo of the system, usually too faint to be observed and/or 
too close to the background quasar. 

Interestingly, our result allows us to strengthen the out- 
flow hypothesis. Considering that 



(i) the mean impact parame ter of Mgll absorbe rs is about 
50 kpc from ai~L* galaxy jZibetti et al.ll2007f ). 

(ii) the mean dust-to-gas ratio of these systems is consis- 
tent with that of L* as a function of redshift, 



this suggests that the gas is originating from the neighbor- 
ing ~ L* galaxy itself. In such a context, outflows appear 
to be the simplest explanation to simultaneously explain 
these two constraints. Dwarf galaxies alone do not satisfy 
the second property. Our result therefore indicates that the 
bulk of Mgll absorbers with W > 1 A and 0.5 < 2 < 1.5 
may trace outflowing gas from galaxies. Another scenario 
would involve low-metallicity gas outflowing or being 
stripped from dwarf galaxies and experiencing a higher 
dust-to-neutral gas ratio due to ionization effects. Such a 
model is however less attractive as it requires ionization 
effects to coincidentally increase the dust-to-neutral gas 
ratio to the value of that of a L* galaxy. 

In some cases, observations of individual galax- 
ies have already pointed out to a link between strong 
Mgll absorbers and outflows: 500 km/s-blueshifted 
Mgll absorption has been observed in the spec- 
tra of 10 out of 14 massive post-starbur st galaxies 
l|Tremonti. Moustakas. fc Diamond-Stank] |2007|). R e cently , 
using integral-field spectroscopy, iBouche et alj (|2007h 
probed H-a emission within 30 kpc of Mgll absorbers 
with Wo > 2 A at z ~ 0.9, and found a 2-cr indication of 
correlation between Mgll rest equivalent width and H-q flux 
based on 14 systems. Using h igh-resolutio n spec troscopy 
of a few strong Mgll systems, iBond et al.1 (|200ll ) showed 
that the shape of certain absorption lines is consistent with 
those expected from galactic winds (but could also be due 
to merging galaxies). While the above analyses dealt with 
individual systems, the statistical results presented in our 
analysis suggest that gas traced by Mgll absorption, with a 
significant range of Mgll rest equivalent widths and over a 
substantial fraction of cosmic time, may often be associated 
with an outflow phenomenon. 
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Figure 4. Relat ion between visual extinction Ay and neutral hydrogen col umn density, N(HI) + 2 N(H2)- D ata points with diamonds 
were taken from ljenkins. Savage, &: Spitzei] lll986h a nd with triangles fromlSnow, Rachford, fe Figoskil j2002h . The line shows the best 
fit obtained from lines of sight within our Galaxy, bv lBohlin. Savage, fc Drake! lll978T ). The data points with error bars show the median 
Nhi and Ay for Mgll absorbers with 1. < Wo < 1.48, 1.48 < Wo < 2.21, 2.21 < Wo < 3.3 by increasing column densities. Blue points 
show direct measurements (with (z) ~ 1) and red points show the extrapolated values to z = taking into account metallicity evolution. 



The diversity of Mgll-selected systems 

Mgll absorption arises in gas spanning several decades 
of neutral hydrogen column density and probes a wide 
range of environments. Locally Mgll absorbers are known 
to trace diverse structures such as the disk of our Galaxy 
Bowen. Blades, fc Pettini 1 19961 ) , in high velocity clouds 



,Bowen, Blades, fc Pettini||l995l; ISayage et all 120001 ) or the 
Large Magellanic cloud (|Weltv et alJll99gf T 



The use of the geometric mean or median hydrogen col- 
umn density and dust-to-gas ratio has allowed us to esti- 
mate number-count weighted quantities rather than Nhi- 
weighted of a Mgll-selected population. Our results aim at 
representing "typical" systems but cannot be extended to all 
objects. The distribution of Nh as a function of Wo seems 
to indicate the existence of two different populations of ob- 
jects: most of the points lie around the dashed lines obtained 
using the geometric mean or median estimates of {Nh ) how- 
ever, others depart by two to three orders of magnitude in 
Nhi from these relations and reach the DLA regime. Such 
systems are likely to have a different nature. First, we can 
observe that the absorbers with the largest hydrogen col- 
umn density do not display the strongest Fell rest equivalent 
width (which better correlates with metal column density 
than Mgll). Furthermore, their implied dust-to-gas ratio is 
necessarily low otherwise the background quasar would be 
heavily extincted. We thus postulate that these systems cor- 
respond, on average, to objects substantially less metal-rich 
such as dwarf galaxies. This statement is in line with the 
mean dust-to-gas ratio and metallic ity estimates of DLAs by 
IVladilo. Prochaska. fc Wolfel (|2008l ). Using about 250 DLAs 



with 2.2 < z < 3.5 detected in SDSS quasar spectra, these 
authors found {A V /N{R1)) ~ 2 to 4 x 10~ 23 mag cm 2 , i.e. a 
dust-to-gas ratio an order of magnitude lower than the mean 
values of Mgll absorbers at z ~ 1. As metallicity evolution 
is n ot expected to be as large between redshifts 1 and 2 
(see Dave fc Opt>enheirnerll2007l or the models summarized 
m iPeroux et alj|2006l ). it indicates that the data points lying 
in the upper part of Fig. [T] have an average dust-to-gas ratio 
significantly lower than those following the median track and 
are likely to trace dwarf-like galaxies. As such, the absorbers 
deviating from the median track and populating the upper 
part of Fig. [T] should correspond to lines-of sight intercept- 
ing galaxies less metal-rich and probably less luminous than 
L* . It will be of interest to test this prediction with available 
or upcoming datasets. A similar conclusion was reached by 
Khare et al. (2007) based on the implied relati on between 
Ay an d N{ZnII) for Mgll selected systems from lYork et al.l 
(2006). Their analysis assumed an SMC dust-to-gas ratio. It 
is worthwile to re-explore their statement for a higher dust- 
to-gas ratio as suggested by the present analysis. 



5 THE Av-Nh RELATION 

As mentioned in the introduction, the 

iBohlin. Savage, fc Drake! (|l978l ) relation, A V /N H ^ 
5.3 x 10 -22 mag cm 2 has been established from measure- 
ments with E(B-V)> 0.1 mag. As our statistical approach 
allows us to be sensitive to reddening values lower by 
an order of magnitude, it is interesting to investigate 
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the validity of the such a relation for (i) lower column 
densities and (ii) in extragalactic environments. In Fig. [4] 
we show measurements o f Ay and Nh for lines- o f-sight 
withi n our Galaxy from Bohlin. Savage, fc Drake! (Il978l ) 
and ISnow. Rachford. fc Figoskil (|2002r i. The blue data 
points show our estimates of (Ay) and the median Nh, 
for three bins of Mgll rest equivalent width. The red 
points show the same quantity but scaled to z = using 
the star formation weigh t ed m etallicity evolution given 
by iDave fc Qppenheimerl (|2007T ) and introduced in the 
previous section. Remarkably, we find that the Bohlin et al. 
relation between Ay and Nh holds down to substantially 
lower hydrogen and dust column densities, regardless of 
ionization corrections. Once evolutionary effects are taken 
into account, the hydrogen and dust content of strong metal 
absorbers at z ~ 1 appear to be in reasonable agreement 
with this relation. 



6 SUMMARY AND OUTLOOK 

The expanded SDSS/HST sample of low-redshift Lyman-a 
absorbers compiled by Rao et al. (2006) has provided us 
with the distribution of hydrogen column densities of strong 
Mgll absorbers. Using this dataset, we have shown that: 

• Due to the nature of the distribution (spanning 3 orders 
of magnitude in Nhi estimates at a fixed Wo), the arithmetic 
mean (Nhi) is effectively sensitive to only a few percents 
of the systems in a Mgll selected sample. Even if (Nhi) or 
Q.hi is computed from a sample of 200 objects, it effectively 
suffers from small number statistics. This fact might be at 
the origin of the discrepancies in SIhi currently debated in 
the literature. 

• Interestingly, despite the large scatter in hydrogen col- 
umn density at a fixed Wo, a well defined relation ex- 
ists between these two quantities: both the median and 
the geometric mean of Nhi show, at the 6-8cr level, that 
(NHi)g.m oc Wq with a ~ 1.8 over a decade in Mgll rest 
equivalent width. 

By combining these hydrogen column density es- 

timations wi th recent reddening measurements from 

Me nard et al.l l|2008l ). we have shown that: 

• while the shape of the mean extinction curve of z ~ 
1 Mgll absorbers is consistent with that of the SMC, the 
inferred dust-to-gas ratio is substantially higher than that 
of the SMC: (A V )/{N(BI)) = 3.0 ± 0.6 x 10~ 22 mag cm 2 . 

• This dust-to-gas ratio does not strongly depend on Wo- 
It varies by less than a factor two for 1 < Wo < 3.3 A 
and suggests that systems in this range have a similar ori- 
gin. Such a property constrasts with the sharp decrease of 
dN/dWo which varies by a factor ~ 30 over the same inter- 
val. This may be due to a transient nature of the absorbing 
gas, such as outflows triggered by star formation. 

• In addition, assuming proportionaly between dust-to- 
gas ratio and metallicity, we have shown that the redshift 
evolution of the dust-to-gas ratio of Mgll systems is in 
agreement with a star formation weigh ted metallicity esti- 
mate (from lrJave fc Oppenheimerll2007l) . Our results there- 
fore confirm the connection between the bulk of strong Mgll 
absorbers and L* galaxies from dust-to-gas ratio considera- 
tions only. 



Interestingly, given that (i) the mean dust-to-gas ratio 
of these systems does not favor LMC/SMC values and 
is consistent with that of L* as a function of redshift, 
and (ii) the mean impact parameter of Mgll absorbers 
is about 50 kpc from a ~ L* galaxy (IZibetti et al.l I2007T ) 
strongly suggest that the absorbing gas originates from the 
nearby ~ L* galaxy. Outflowing gas appears as the most 
simple explanation regarding the nature of the majority 
of strong Mgll absorbers. An alternative scenario would 
involve low-metallicity gas originating from dwarf galaxies 
and experiencing a higher dust-to-neutral gas ratio due to 
ionization effects. Such a model is however less attractive as 
it requires ionization effects to coincident ally increase the 
dust-to-neutral gas ratio to the value of that of a L* galaxy. 

Our results made use of geometric mean and/or median 
estimates. They are aimed at representing the bulk of Mgll 
absorbers and are not sensitive to outliers in the distribu- 
tions of hydrogen column densities. We have shown that the 
systems populating the upper part of Fig [T] and correspond- 
ing to DLAs are necessarily less dusty and likely less metal- 
rich than the bulk of Mgll systems indicated by the median 
track. This illustrates that, even if an overlap exists between 
the two populations, the majority of DLA-selected or Mgll- 
selected systems do not probe the same type of structures. 

It will be of particular interest to repeat such an analy- 
sis for metallicity estimates of absorber systems. In addition, 
accessing several transistions would allow us to better con- 
strain the ionization fraction. 

Finally, we have shown that our measurements allow 
us to probe the relation between Ay and Nh, i-e. between 
dust and hydrogen column densities, in regimes that were 
prev iously unexplored. In par ti cular, we have shown that 
the iBohlin. Savage, fc Drake] (|l978l ) relation measured 
from lines-of-sight in our Galaxy is also satisfied by Mgll 
absorbers at z ~ 1, for column densities and reddening 
values significantly lower than previously probed. As it 
relates the amount of dust to the neutral hydrogen column 
density, this statement is valid regardless of the level of 
ionization. 

If the bulk of Mgll absorbers are tracers of outflowing 
gas, they may provide us with an interesting view of delayed 
star formation around galaxies. Assuming a wind velocity of 
about 300 km/s, i.e. about 300 kpc/Gyr, the presence of a 
strong Mgll absorber is then typically related to a burst of 
star formation which occured about 150 Myr ago. The ob- 
served properties of Mgll absorption lines (incidence, red- 
shift evolution, velocity width, number of components) can 
provide us with observational constraints on star formation 
in the Universe, up to high redshifts where emission studies 
are more difficult. 
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